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Free-standing exfoliated monolayer graphene is an ultra-thin flexible membrane, which exhibits 
out of plane deformation or corrugation. In this paper, a technique is described to measure the 
band structure of such free-standing graphene by angle-resolved photoemission. Our results show 
that photoelectron coherence is limited by the crystal corrugation. However, by combining surface 
morphology measurements of the graphene roughness with angle-resolved photoemission, energy- 
dependent quasiparticle lifetime and bandstructure measurements can be extracted. Our measure- 
ments rely on our development of an analytical formulation for relating the crystal corrugation to 
the photoemission linewidth. Our ARPES measurements show that, despite significant deviation 
from planarity of the crystal, the electronic structure of exfoliated suspended graphene is nearly 
that of ideal, undoped graphene; we measure the Dirac point to be within 25 meV of Ef- Further, 
we show that suspended graphene behaves as a marginal Fermi-liquid, with a quasiparticle lifetime 
which scales as (E — Ef) -1 , comparison with other graphene and graphite data is discussed. 



PACS numbers: 73.22.Pr,68.65.Pq 
I. INTRODUCTION 

The recent availability of monolayer-thick two- 
dimensional crystals such as graphene, BN, and BSCCO 
has generated widespread interest in the physics and ma- 
terials science communities. In the case of graphene, in 
particular, the two dimensional nature of the crystal in 
combination with its unusual massless Dirac fermions de- 
termines a host of intriguing and unique transport phe- 
nomena, including graphene's half-integer quantum Hall 
effect (HE) and non-zero Berry's phaseP^ Unlike most 
metals, undoped graphene has a Fermi surface which con- 
sists of a set of 2 inequivalent points in momentum- space. 
Thus, at zero temperature and zero doping, the density 
of states at the Fermi level vanishes. In combination with 
the linear dispersion of low energy charge carriers, this 
vanishing density of states is expected to lead to unusual 
band-renormalization effects that are not seen in Fermi- 
liquid systems such as unusually high electron-electron 
coupling. Motivated by interest in these unusual prop- 
erties, several theoretical and experimental studies have 
investigated the electronic properties of grapheneP 

Angle-resolved photoemission spectroscopy (ARPES) 
is the experimental method that is most frequently used 
to probe the electronic structure of crystals. However, 
so far, the majority of ARPES studies of graphene have 
been conducted on epitaxial graphene, which has been 
grown on a variety of substrates such as SiC, Ru, Ni 
and IrPK^ Epitaxial graphene is ideal for photoemission 



experiments, but, due to the interaction between the epi- 
taxial graphene monolayer and the substrate, the band 
structure is often distorted such that the Dirac point 
shifts away from the Fermi energy, thus changing the 
quasiparticle dynamics. In an effort to minimize the ef- 
fect of substrate interaction on epitaxial graphene, recent 
ARPES studies have focused on several multilayer sys- 
tems, such as intercalated graphite^ and graphene grown 
on the C face of SiC These layered systems consist of 
multiple stacked graphene sheets that are substantially 
electrically isolated, thus resulting in an electronic band 
structure that mimics that of suspended exfoliated single- 
layer graphene. However, despite its scientific and tech- 
nological importance, exfoliated graphene has been the 
subject of only a limited number of ARPES studiesj^^ 
despite the fact that it remains the best choice for de- 
vice physics, as it is easily backgated and has the highest 
measured mobility.^ 

Several obstacles impede measurement of the band- 
structure of exfoliated graphene. One difficultly arises 
from the fact that available single-layer exfoliated 
graphene flakes are typically less than 20 jam in size, 
thus precluding the use of standard ARPES systems, 
which require samples to be several mm in size. Hence, 
most information regarding low-energy occupied states 
in exfoliated graphene has been obtained indirectly 
from electrical-transport measu rements^ or directly by 
optical-probing techniques! 16 * 17 * These techniques exam- 
ine the bandstructure generally within leV of the Dirac 
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point and do not directly provide momentum resolution. 
For photoemission the limitation in size can be overcome 
by working with high lateral-spatial-resolution prob es 
such as those available using spectromicroscopyP^^ A 
second major impediment to photoemission studies is due 
to the fact that graphene is an ultrathin crystal. This 
ultrathin property has, in turn, two important conse- 
quences for photoemission studies. The first is the trans- 
parency of monolayer graphene to UV photons and pho- 
toemitted electrons, which causes a strong background 
photoemission signal if the monolayer graphene is in close 
physical proximity with a substrate. 14 The second is that 
exfoliated graphene is not atomically flat, but is known to 
deform locally, a result shown through AFM, STM, elec- 
tron microscopy, and electron scattering results ! 19 " 23 l It 
has been argued that the deformation is due to the fact 
that monolayer-thick graphene has soft flexural modes 
leading to ready bending of the graphene. The pres- 
ence of a supporting substrate or scaffold can, to a cer- 
tain degree, stabilize height fluctuations in the graphene 
layer, but corrugations in the underlying supporting sub- 
strate are transferred in part to the graphene due to the 
reduced stiffness of this material. Additionally, intrin- 
sic corrugations that cannot be attributed to interaction 
with the substrate were recently observed in supported 
graphene Further, in a recent low energy diffraction 
study, we demonstrated that even graphene suspended 
over etched cavities exhibits corrugation, which appeared 
to have been intrinsic in originPH 

Thus, in general, two dimensional crystals produced by 
exfoliation may show significant local curvature, mani- 
fested as corrugation and ripples. This corrugation is 
known to affect not only the electronic and transport 
properties of the material, but can also have a major 
impact on photoemission results. In particular, the the- 
ory of ARPES was developed for single-crystal atomi- 
cally flat surfaces and relies on the fact that momentum 
perpendicular to the surface is conserved in the photoe- 
mission process. On such perfectly ordered crystals the 
photoemission lineshape is directly related to the spectral 
function of the electronic state being probed, from which 
information about many-body physics can be extracted. 
The corrugation in thin sheets of layered materials breaks 
this symmetry and obscures the intrinsic many-body ef- 
fects. 

In this paper, we present a systematic approach 
to account for such corrugation-induced broadening in 
ARPES on thin films. By combining our photoe- 
mission results with detailed information about sur- 
face morphology obtained from prior electron-microscopy 
measurements 21 taken in-situ on the same samples we are 
able to quantify the influence of corrugation on spectral 
broadening. We go on to describe a method to discount 
the effect of surface corrugation from ARPES measure- 
ments to reveal the intrinsic many-body physics present 
in graphene. Our results show that suspended graphene 
behaves as a marginal Fermi-liquid with an anomalous 
quasiparticle lifetime which scales as [E — Ep)~ x . 



II. EXPERIMENT 

Our measurements used the Spectroscopic Photoemis- 
sion and Low Energy Electron Microscope (SPELEEM) 
at the Nanospectroscopy beamline at the Elettra Syn- 
chrotron light source. 24 The SPELEEM is a versatile 
multi-technique microscope that combines low energy 
electron microscopy (LEEM) with energy-filtered X-ray 
photoemission electron microscopy (XPEEM). The mi- 
croscope images surfaces, interfaces and ultra-thin films 
using a range of complementary analytical character- 
ization m ethod s, which have been described in detail 
previously! 25 ! 26 * When operated as a LEEM, the micro- 
scope probes the specimen using elastically backscattered 
electrons. LEEM enables high sensitivity to surface crys- 
talline structure and, due to the favorable backscattering 
cross-sections of most materials at low energies, allows 
image acquisition to be obtained at video frame rate. The 
lateral resolution of the microscope for LEEM imaging is 
currently below 10 nm. In XPEEM mode, the specimen 
is probed using the beamline photons, provided by an un- 
dulator source; thus, the technique is sensitive to the lo- 
cal chemical and electronic structures. Laterally resolved 
versions of synchrotron based absorption (XAS) and pho- 
toemission spectroscopy (XPS) are possible. The lateral 
resolution in XPEEM approaches a few tens of nmP 

Along with real-space imaging, the SPELEEM mi- 
croscope is capable of micro-probe diffraction imaging, 
i.e. laterally restricted low energy electron diffraction 
(/i— LEED) and angle resolved photoemission electron 
spectroscopy (fi— ARPES) measurements when probing 
with electrons and photons, respectively. In diffraction 
operation the microscope images and magnifies the back 
focal plane of the objective lens. In ARPES mode, the 
full angular emission pattern can be imaged on the de- 
tector up to a parallel momentum of ~ 2A ; at larger 
parallel momentum the transmission of the microscope 
decreases. All diffraction measurements are restricted 
to areas of ~ 2 fim in diameter, which are selected by 
inserting a field limiting aperture into the first image 
plane along the imaging-optics column of the instrument. 
Thus, the microscope enables measurements on samples 
that are homogeneous over areas of a few square microns. 
The energy resolution of the SPELEEM in diffraction 
imaging, such as ARPES and photoelectron-diffraction 
measurements, is 300 meV and the transfer width of the 
microscope when operated in LEED mode is 10 nm. 

Graphene samples were extracted by micro-mechanical 
cleavage from Kish graphite crystals (Toshiba Ceramics, 
Inc.) and placed onto an Si02-thin-film layer on an Si 
substrate, which was previously patterned with cylindri- 
cal cavities to a depth of 300 nm, as described in Ref. [21] 
The use of planar processing of this substrate allowed us 
to suspend areas of the graphene films without the use 
of further photolithographic techniques, which would in- 
troduce contaminants to the graphene sheets. Graphene 
samples with lateral sizes from 10 to 50 /im were placed 
in contact with Au grounding strips deposited on the sur- 
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FIG. 1. (Color Online) (Left) Schematic drawing of our 
suspended-graphene sample configuration. (Right) Opti- 
cal micrograph of sample containing suspended monolayer 
graphene (MLG) and few- layer graphene (FLG). (Bottom) 
Artists rendering of corrugated graphene crystal (height fluc- 
tuations not to scale). 
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FIG. 2. (Color online) Raw ARPES data along symmetry di- 
rections in Brillouin Zone for graphene and graphite. (a)Si02 
supported graphene (/kj=90 eV). (b) Suspended graphene 
(ftcj=84 eV). (c)Kish graphite (/kj=90 eV). (d) Suspended 
graphene (hou=50 eV). Inset shows 2D graphene Brillouin 
zone. 



face via thermal evaporation through a shadow mask. A 
sketch of the sample configuration is shown in Fig. 1 
along with an optical micrograph. 

The SPELEEM instrument used to collect data has the 
important advantage of having a sufficiently high spatial 
resolution to guarantee that we are measuring a single 
crystal sample of monolayer graphene and that all of 
the measured spectral intensity is derived from a fully 
suspended region. This capability is necessary since the 
suspended regions are approximately 5 jam in diameter 
and, therefore, cannot be resolved with conventional pho- 
toemission instruments, which employ spatial averaging 
techniques that collect data over surface areas of sev- 
eral square millimeters. The potential to combine both 
photoemission and electron scattering measurements is 
essential for our experiment since it allows us to measure 
bandstructure and surface morphology on the same sam- 
ples. We note, additionally, that a similar instrument was 
recently used in a study, which examined the morphology 
and electronic structure of epitaxial graphene grown on 
PtP 

After preparation the samples were placed into a UHV 
chamber with a base pressure of 2 x 10 -10 mbar, and 
the surface cleaned via low energy electron irradiation to 
eliminate adventitious hydrocarbon molecules adsorbed 
during prior atmospheric exposure. 21 All graphene sam- 
ples were characterized with LEEM before investigation 
with ARPES and LEED. For each sample, LEEM was 
used to locate sample areas of interest and to determine 
film thickness with atomic resolution by mea suring in- 
tensity modulations in the LEEM I-V spectraP- 1 | 28 | 2Q | 

ARPES data at multiple photon energies were ob- 



tained on the suspended areas of the graphene film. Only 
regions of uniform thickness were considered. In order to 
elucidate the role of surface corrugation and substrate 
influence, comparative experiments were also carried out 
on corresponding regions where the film was supported 
by the Si02 substrate. This surface has been recently 
carefully calib rated by prior STM and electron-scattering 
measurements 112120122] j n addition, ARPES measurements 
were made on Kish-graphite flakes that were present on 
the same substrates. As graphite is a well understood and 
commonly studied system, these measurements provided 
a useful point of comparison for our graphene measure- 
ments. Photoemission from graphite is, in some respects, 
similar to that from graphene because of the stacked- 
layer nature of the former. However, the physics near the 
Dirac point is significantly different owing to the fact that 
the multilayer stacking in graphite breaks the symmetry 
between A and B sublattices, which results in two dis- 
persing branches, such that low energy excitations do not 
have the simple linear dispersion relation that is found 
for graphene. 



III. RESULTS 

Photoemission spectra were measured from two sam- 
ples with differing degrees of surface corrugation and sub- 
strate interaction, that is, on suspended and substrate- 
supported graphene. Previous LEED measurements have 
shown that the horizontal correlation length increases 
from 24 nm to 30 nm in measurements taken on sup- 



4 



ported and suspended samples, respectively.^ In addi- 
tion, ARPES data were collected at room temperature 
over the entire surface Brillouin zone (SBZ) from 0.5 
eV to -8 eV (energy referenced to Ep), for monolayer 
graphene and graphite, using a range of photon ener- 
gies. Figure [2] shows ARPES spectra taken from a sam- 
ple supported by and in contact with the Si02 surface 
and a sample that was suspended over the 5 jira wells 
shown in Fig. [I] For comparison, the raw ARPES spec- 
trum from Kish graphite is shown as well. The data 
show dispersion along 3 symmetry lines in the SBZ. As 
expected from the reduced corrugation, as well as the 
absence of any substrate interaction, the ARPES data 
for suspended graphene show a dramatic improvement in 
quality as compared to the data for supported graphene. 
Additionally, there is a very broad, parabolically dis- 
persing peak centered at the T point at a binding en- 
ergy of ~8 eV in the data taken on supported graphene. 
This feature has been previously attributed to photoe- 
mission from the amorphous Si02 substrate^ and is not 
observed in the spectrum taken on suspended graphene. 
Although the substrate is only 300 nm below the sus- 
pended graphene, any background electrons emmited at 
this height will be significantly defocused in the electron 
optics of SPELEEM microscope. Additionally, due to 
the grazing incidence angle of the photon beam (16°), 
the bottom of the cavity is not fully illuminated as the 
cavity edge casts a shadow, which further reduces the 
photoemission signal from the substrate. 

In the vicinity of the K points, a conical dispersion 
is observed centered at the K point on the suspended 
graphene spectrum. At ~ 1 eV below the Fermi level 
a trigonal-warping deviation from angular isotropy be- 
comes clearly noticeable. Measurements taken through 
the K point and in the direction parallel to the TM direc- 
tion (vertical direction) show two symmetric dispersing 
branches forming the two sides of the Dirac cone. The 
band structure can be made significantly sharper (see 
Fig. [3| by taking the second derivative along each mo- 
mentum direction. In this case, use of the second deriva- 
tive allows easier determination of the Dirac point with 
respect to the Fermi level. Figure shows the lin- 
ear best fit to the two branches as well as the location 
of the Fermi level. From the fit, we find that the Dirac 
point is within 25 meV of E F (E D = — 9 ± 25meV). 
Thus, the sample is minimally doped due to the prepara- 
tion procedure used here, which did not involve any pho- 
tolithographic or chemical-transfer techniques. In con- 
trast, the Dirac point previously measured by our group 
on a supported sample was found to be ~300 meV be- 
low the Fermi level, which was attributed to doping by 
interaction with charged impurities in the Si02 layer 

For comparison with results on a known photoemis- 
sion materials system, graphite spectra were taken at two 
photon energies (86 and 76 eV) along the same (vertical) 
direction through the K point; these results are shown in 
Fig. [4] The dispersion obtained at huu=86 eV is clearly 
symmetric about the K point. At this photon energy we 
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FIG. 3. (Color online) (a) ARPES intensity through K point 
along TM (k y ) direction in suspended monolayer graphene 
(huo = 50eV). (b) Smoothed second derivative image of dis- 
persion shown in (a), (c) Smoothed second derivative ARPES 
intensity through K point along TK (k x ) direction, (d) Ex- 
tracted dispersion from (c) Inset shows graphene Brillouin 
zone. Red solid (dashed) line indicates k y (k x ) direction 
through K point. 
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FIG. 4. (Color online) Dispersion along the vertical (k y ) 
direction through the K point for graphite obtained at pho- 
ton energies of (a) huo = 86 eV and (b) huo = 76 eV. (c-d) 
Smoothed second derivative images of spectra shown in (a) 
and (b), respectively. Inset shows graphite surface Brillouin 
zone. Solid red line indicates k y direction through K point. 
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FIG. 5. (Color online) ARPES intensity along TK direction 
in suspended monolayer graphene. Graphene photoemission 
taken with photon energies of (a) hw = 50eV and (b) hw = 
84eV. (c) HWHM of MDCs as a function of binding energy 
taken from (a) and (b). Inset shows sample MDCs taken 0.7 
eV below Ef as indicated by dashed red lines in (a) and (b). 



can resolve the splitting of the tt state into bonding and 
antibonding bands, with the two bands separated by ~ 
0.12 A . The bands themselves are approximately 0.1 
A in width. In the spectrum taken at 76 eV the two 
peaks are nearly degenerate. Again, the second deriva- 
tive allows for easier determination of peak locations. 

Figure [5] shows the graphene dispersion taken along the 
TK direction through the K point. Comparative mea- 
surements were made at two different photon energies 
(hw = 50, hw — 84) and are shown in Fig. [5|a) and (b), 
respectively. In this direction, only one branch of the 
dispersion can be seen as the difference in phase between 
electron waves emitted from the A and B sub-lattice sites 
results in complete destructive interference! 3 ^ Thus, this 
is a convenient direction along which to measure pre- 
cisely the dispersion in the vicinity of the Dirac cone. 
The inset to Figure [5|c) compares momentum distribu- 
tion curves (MDCs) taken at a binding energy of 0.7 eV 
for the suspended-graphene spectra at both photon en- 
ergies. The data in Fig. [5|c) show that the width of the 
huj=S4: eV MDC is significantly larger than the hw=50 
eV MDC (0.17 A" 1 vs 0.12 A" 1 ). 

Additionally, there is a slight asymmetry in all three 
MDCs as additional spectral weight is present on the 
right side of the peak (at higher values of k x ). The 
background signal decreases and the peaks become nar- 



rower for 50 eV photons as compared to 84 eV photons. 
Specifically, in the 0-4 eV range (referenced to Ef), the 
MDC width increases monotonically from 0.1 to 0.2 A 
and from 0.15 to 0.3 A for data collected with 50 eV 
photons and 84 eV photons, respectively. In contrast, 
MDCs taken along the same direction (rK) on supported 
graphene are significantly broader 14 and show almost no 
dependence on binding energy; they are -0.5 A m 
width from the Fermi energy to -4 eV binding energy. 
Thus, spectral features are sharpest for suspended sam- 
ples measured with lower photon energy. 



IV. DISCUSSION 
A. Comparison of Graphite and Graphene Results 

As shown in Fig. [5j variation in photon energy results 
in changes to the linewidth of the graphene photoemis- 
sion spectra, which can be exploited to sharpen the spec- 
trum. In explaining these results on graphene, it is useful 
first to examine the effect of photon energy variation for 
the case of graphite. The differences in the measured 
photemission spectra of graphite taken at hw=76 eV and 
hw=86 eV shown in Figs. [4^a) and (b), respectively, 
are easily understood by considering the 3 dimensional 
band structure of graphite. In particular, according to 
the standard model of photoemission, variation in hw al- 
lows one to access a range of initial states with different 
k z W^ Using the free-electron approximation for t he fi- 
nal state allows calculation of k z of the initial stateP^^ 
Thus, in the case of the data shown in Fig. [4| photoe- 
mission obtained at hw = 86 eV corresponds to k z = 
(the TKM plane), while hw = 76 eV accesses k z = 0.3c* 
which is nearer the AHL plane. Since k z is changed, the 
clear double band feature seen for k z = changes as the 
graphite band structure varies alo ng k z in accord with 
the known graphite band structure! 30 * 32 " 34 ! 

Consider now the effect of changing photon energies 
for the case of graphene photoemission. Since graphene 
is truly a 2D crystal, the initial states in the valence band 
are highly localized along the z direction. Thus, the Bril- 
louin zone is strictly 2 dimensional and the electronic 
strucure is essentially k z independent. Comparison with 
photoemission from surface states is useful, since they are 
also localized in 2D. 35 However, the role of evanescent de- 
cay into the bulk, which is important for surface states 
in metals and results in a partial k z dependence^ such 
as surface resonance, is absent in graphene and, thus, 
we may treat the initial state as independent of photon 
energy. In fact, as seen in Fig. [5j changing the photon 
energy in the case of graphene causes only a change in 
the overall linewidth and does not affect the measured 
bandstructure. As will be discussed below, the differ- 
ence in the width of ARPES features between spectra 
obtained at hw =50 and hw =84 is a consequence of the 
surface roughness of the graphene samples. Since elec- 
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trons in graphene propagate on a locally curved surface, 
the usual momentum conservation rules in ARPES must 
be modified and a photon-energy-dependent broadening 
term is introduced. 



B. General Considerations 

In standard many-body ARPES theory, the intensity 
of the photoemission signal is proportional to the spectral 
function, A(k,cj): 



A(k,£j) 



Im[E(k,c 



(uj - cj k - Re[E(k, uj)}) 2 + Im[(E(k, uj)} 2 



(1) 



where uj = E — Ep and k are binding energy and mo- 
mentum, respectively, and uj^ is the single-particle dis- 
persion. The real and imaginary parts of the self-energy, 
E(k, cj), represent renormalization of the bare-bands and 
scattering rate, respectively. To obtain the full expression 
for the photocurrent, the above function is then multi- 
plied by energy and momentum-preserving delta func- 
tions, 5(ki — k/ — G)S(Ei — Ef — W), where G is a recip- 
rocal lattice vector and i and f label the initial and final 
states, respectively, and W is the work function of the 
material. 

However, one major complication to this approach 
arises for the case of suspended graphene since the mo- 
mentum preserving function, S(ki — kj — G), is only a 
precise delta- function if the system under investigation 
is atomically fiat. While this is the case for the majority 
of single-crystal samples probed with ARPES, including 
the Kish graphite described above, exfoliated monolayer 
graphene, as is discussed in the Introduction, has sig- 
nificant deviations from planarity, ranging from 1 to 10 
j^pO This corrugation introduces an additional broaden- 
ing mechanism into the ARPES spectrum, which can be 
as large as, or larger than, the intrinsic broadening repre- 
sented by Im[E(k, lj)]. Thus, in order to extract the true 
self-energy of carriers in the crystal, such corrugation- 
induced broadening must be taken into account. The 
MDCs are best fit by a convolution of A(k, uj) with a func- 
tion that represents broadening due to surface roughness. 
Thus, as will be shown below, at fixed cj, photoemission 
intensity as a function of kn can be expressed as: 



7(k||)oc J d\{S k A\\)A{\\-^u) 



(2) 



where ky = k^y — kjy and S^ ± represents corrugation- 
induced broadening. S\^ ± , the surface structure factor, is 
a function of the surface geometry of the sample and is 
also generally dependent on the change in perpendicular 
momentum from initial to final state, k± = ki± — kf±. 
We note that several prior studies have examined t he ef- 
fect of surface roughness on ARPES measurements.^^ 
In these prior studies, the roughness considered was due 



to discrete height variations caused by monatomic steps, 
rather than the continuous undulations of a thin film. 
Thus, the broadening in spectral features measured by 
ARPES was attributed to increased electron scattering 
rather than a variation in the phase of photoemitted elec- 
trons induced by local height fluctuations. In our exper- 
iments on suspended graphene samples, the surface mor- 
phology is carefully measured simultaneously with the 
ARPES measurements presented here, thus allowing us 
to determine S\^ ± independently.^ 

Finally, note that the surface corrugation of the 
graphene sheets will also alter the bandstructure by in- 
ducing a change in the local potential proportional to the 
square of the curvature. Thus, the ripples act as scatter- 
ing centers, which will decrease lifetime and potentially 
change the Fermi velocity. These effects are contained 
in A(k, uj) and will also be present in the ARPES data. 
However, such effects are distinct from that described by 
Sk ± , which represents decoherence as electrons pass from 
a curved 2D space to free space. 



C. Corrugation Broadening 

Corrugation broadening can be treated by consider- 
ing the equation that describes photoemission from a 
Bloch state in the graphene sheet into a free-electron 
state above the crystal. Using the standard tight-binding 
approach to describe the initial state: 



W = ^E^ R E CSWr-R-r,-) (3) 

V7V R j=A,B 



l/jl 



we obtain the following matrix element for excitation into 
a free electron final state: 



Moc(ki-A) ^ ki e^ kf - r ^e i(ki " kf),R 0(k f ) (4) 

j=A,B R 

where ki is the initial pseudo-momentum of a valence- 
band electron and kf is the final-state momentum (for a 
full description and definitions of symbols see Appendix) . 
Equation [3] describes an initial state with precise mo- 
mentum at a fixed binding energy. For an atomically- 
flat crystalline 2D surface the position vectors can be 
expressed as R = niai + 712^2, where the ni are inte- 
gers and the primitive lattice vectors in the xy plane. 
In this case, the sum over R in Eq. [5J ^2 e *(ki-k f )-R^ 
is zero unless k^ — kfy = G, where G is a recip- 
rocal lattice vector. This condition is, thus, a state- 
ment of the momentum conservation discussed above, 
S(ki — kf — G). If, however, z is allowed to vary con- 
tinuously as a function of position along the surface, so 
that R = nidLi + 7128-2 + Ax + Ay + J^, with z no longer 
constant, the summation in Eq. [4] is not as readily cal- 
culated. Perfect phase cancellation away from reciprocal 
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lattice vectors does not occur, resulting in non-zero pho- 
toemisison intensity when k^ — kf y ^ G. 

Summations such as the one in Eq. [4] are encoun- 
tered in the theory of LEED on rough 
fact, in many respects the formal analysis of LEED re- 
sults bears many similarities to that of ARPES. In a 
prior study using one-photon photoemission and high- 
resolution LEED applied simultaneously to surface states 
of Cu(100) and Cu(lll), it was demonstrated experimen- 
tally that the photoemission linewidth and the width of 
the LEED-spot profile are correlated linearly.^ In par- 
ticular, for LEED one measures the diffraction structure 
factor, S(k) oc | e lk R | 2 where, as in the case of photoe- 
mission, k is the total momentum transfer, k = ki — kf , 
and the sum is over atomic positions, R, on a surface. 
In addition, for ARPES transition probability is propor- 
tional to the square of the matrix element; thus, the same 
structure factor, 5(k), is applicable. Thus, LEED theory 
can guide our analysis. 

The structure factor, 5'(k), can be calculated with in- 
formation about the average properties of the surface, de- 
scribed by three variables: horizontal correlation-length, 
£, RMS height variation, w, and a dimensionless parame- 
ter, a, termed the "roughness exponent," which describes 
surface roughness on length scales smaller than All 
three parameters can be extracted from real-space infor- 
mation about the surface by computing the height-height 
correlation function, which is used in a variety of thin film 
measurements, including those on graphene and other 
surfaces, and is defined as H(r) = (\z(ro + r) — z(ro)\ 2 ). 
As is shown in the appendix, 5(k) = 5(fcj_,ky) is in- 
timately related to H(r) as the Fourier transform of 
e -^k±H(r) _ Thus, the average parameters that character- 
ize a given rough surface (w, £, and a) and determine the 
form of H(r) also determine S(k± : ky). Hence, with these 
parameters, it is possible to compute the summation in 
Eq. [4j In fact, previously reported measurements us- 
ing low-energy electron microscopy and low-energy elec- 
tron diffraction have determined these parameters to be 
a = 0.54 ± 0.02, w = 1.99 ± 0.15 A, and f = 30 ± 0.3 
nm for the same suspended graphene samples used in 
this study!^ Although the functional form of S f (fej_,ki|) 
is complex, the width of S^^ky) (i.e. for k± fixed) in ky 
space has a simple dependence on k± and the parame- 
ters describing the surface roughness. In particular, the 
width, is proportional to (k^w) 1 ^ /£, which explains 
the decrease in experimental linewidth with decreasing 
k± shown in Fig. [5] For fitting purposes it is useful 
to have the exact functional form of 5^(ky). Yang, et 
al have shown that for (wk±) 2 ^> 1 the form is purely 
diffusive and can be expressed asP^ 

S^(ky) = ^/(wk^-)F a (^/(wk^-) 

F a (Y) = [ XdXexp(-X 2a )J (XY) (5) 




-0.4 -0.3 -0.2 -0.1 ' -0.2 -0.1 0.1 0.2 



k x k y 

FIG. 6. (Color online) (a) Graphene bandstructure over first 
BZ. (b) Contour lines drawn along constant binding energies 
in the vicinity of the K point (binding energies indicated in 
eV). (c-d) Dispersion along k x and k y , respectively (along 
dashed blue lines in (b)). 



D. Intrinsic Broadening 

It is straightforward to introduce intrinsic initial-state 
broadening into our ARPES description by replacing our 
initial state wavefunction, ip^ with a sum over multiple 
momentum states, J^ak^k? where the a ki are complex 
coefficients related to the spectral function by |a k .| 2 = 
A(k, uj). Our transition matrix then becomes a sum, 
M = J^akM k , over multiple matrix elements weighted 
by the complex coefficients <2k, where the M k are the 
original transition matrix elements defined in Eq. [4] 
Again, using Fermi's golden rule we find that the tran- 
sition probability is proportional to the square of this 
sum. 



Ioc\M\ 2 = |^a fc M k | 2 = ^|a k | 2 |M /c | 2 +^ a k a k ,M k *M k/ 

k k k^k' 

(6) 

As shown in the appendix the k ^ k 7 sum can be safely 
neglected due to random phase cancellation and we arrive 
at the final expression for the full photoemission intensity 
expressed in Eq. [2] 

Finally, we note that, in general, the linewidth (r m ) 
measured in ARPES from well prepared, atomically fiat 
surfaces is a function of the initial state or photohole 
linewidth (1^) as well as the linewidth of the final state 
or photoelectron (Tf). However, for the case of 2D states 
such as surface states in metals or thin films such as 
graphene, there is no dispersion with k± and T m = Ti^i 
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E. Analysis of Spectra and Discussion 

Figure [3] shows a plot of the dispersion obtained along 
the TK direction in the vicinity of the Dirac point. The 
average Fermi velocity, derived from the slope of uj vs 
k x is 1.07 ± 0.05 x 10 6 m/s. This value is in excellent 
agreement with results obtained by IR m easur ements on 
undoped supported exfoliated graphene ! 44 1 45 1 Addition- 
ally, the dispersion along TK is linear with no deviations 
from linearity within our experimental uncertainty. As 
discussed above, despite the roughness induced broaden- 
ing in the spectrum, the dispersion curve is easily ex- 
tracted from the raw ARPES data by taking the second 
derivative of the ARPES intensity along the momentum 
direction. However, determining the intrinsic width of 
spectral features requires a deeper analysis. 

Our prior measurements of the surface corrugation in 
suspended graphene allow us to extract the intrinsic elec- 
tronic structure from our ARPES data. The procedure 
for this fitting is as follows: first, Sk ± is determined 
from our surface morphology measurements and used as 
a constant parameter, then A(kn) is varied until the con- 
volved function, /(ky), represents a good fit to the ex- 
perimental data. Although a full deconvolution is, in 
principle, possible it is much more straightforward to be- 
gin with an assumption for the functional form of A(kn) 
and systematically vary the parameters until a good fit 
is found. The functional form of A(ky) is assumed to 
be a Lorentzian, the most commonly used photoemission 
lineshape, which results from the k-independent approx- 
imation for Im[E(k, a;)] P 

In carrying out this procedure, we introduce two ad- 
ditional simplifications. First, note that we are examin- 
ing the region in k-space within the first Brillouin zone 
along the TK (k x ) direction in the vicinity of the K point. 
As shown in Fig. [6| the tt state disperses rapidly along 
k x in this region, but relatively slowly along k y since 
doj/dk y =0 at k y =0. Thus, although Eq. [6] describes a 
2D convolution, it is possible to replace the required 2D 
k x k y integral with a ID integral along k x . Second, we 
note that most of the MDCs considered here have an 
asymmetric peak shape, with additional spectral weight 
in the uj > vp\k\ region of the curve. Possible reasons for 
this asymmetry are discussed in a separate paragraph be- 
low. For fitting purposes, this additional spectral weight 
was not considered and the best fit was obtained by im- 
posing a momentum cutoff within 0.1 A of the peak 
position on the uj > vp|k| side of the curve. Figure [7^a) 
shows a representative curve from the huj=50 eV data 
taken 0.7 eV below the Fermi level along with a best 
fit. Note that the lineshape of this curve provides an ex- 
cellent fit to the experimental data. Figure [T^b) shows 
the two independent contributions to the linewidth: the 
corrugation-induced broadening and the intrinsic broad- 
ening. In order to cross-check that the convolution proce- 
dure accurately captures the photon-energy dependence 
of the photoemission process, the same fitting procedure 
was repeated on data obtained with a photon energy of 




FIG. 7. (Color online) Intrinsic width of ARPES features for 
Suspended monolayer graphene. (a) Example of MDC fitting. 

(b) Two independent contributions to broadening. Red line: 
intrinsic linewidth of initial state. Blue line: broadening due 
to corrugation at /kj=50 eV (solid) and /kj=84 eV (dashed). 

(c) Inverse lifetime as a function of binding energy for Huj=50 
(blue) and /kj=84 eV (red). Inset shows best fit line to in- 
trinsic width vs. binding energy for hou=50 data in vicinity of 
Fermi level. 

7kj=84. At this photon energy, k± = 4.27A 1 and, ac- 
cording to Eq. |6| the width of Sk ± is nearly twice as 
large as it is at Hlj = 50 eV. However, as expected, the 
intrinsic linewidth extracted from the fitting procedure 
is the same for data obtained with both photon energies. 
A comparison of the self-energy extracted from the two 
data sets is shown in Fig. [7^c); the two resulting curves 
are the same, within experimental error, thus confirm- 
ing the photon-energy dependence given in Eq. [6] and 
lending further support to our approach. 

To make our observations quantitative and enable com- 
parison with other work, we perform a linear fit of the 
intrinsic width versus binding energy, I\ = C6+/3(E—Ef). 

From this fit, we find a = 0.002 ± 0.005 A" 1 and 
P = 0.039 ± 0.01 A^eV- 1 . As expected, the value of 
a is within experimental uncertainty of zero since ex- 
cited states just above the Fermi level should be very 
long lived. Consider now the parameter f3 that de- 
scribes the increase in inverse quasiparticle lifetime with 
increasing binding energy. The lifetime is related to I\ 
by r = l/(2r^ j p). Thus, our measured value can be 
reexpressed as j3= 0.78±0.02 fs _1 eV -1 , so as to en- 
able ready comparison with prior measurements of the 
same quantity on graphite and exfoliated graphene. For 
graphite, /3 has been measured by femtosecond pho- 
toemission to be an order of magnitude smaller, viz. 
0.029 fs _1 eV -1 } 4 ^ while STS measurements of exfoli- 
ated graphene on graphite have produced an intermediate 
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value of (f3 = 0.11 fs _1 eV _1 )!^A reasonable explanation 
for this discrepancy is the greater out-of-plane corruga- 
tion of suspended graphene, which has been predicted 
to be the largest contribution to electron scattering in 
rough graphene sheet iP^M Indeed, such roughness con- 
stitutes short-range correlated disorder, which has also 
been shown theoretically to lead to scattering rates which 
scale linearly with uj in graphene!^ 

Comparison can also be made with results obtained 
on epitaxial graphene grown on SiC. In such a system 
the Dirac point is ~ 0.5 eV below the Fermi level which 
changes the quasiparticle dynamics resulting in a non- 
linear behavior for I\ vs uj. In particular, it has been 
shown that electron-plasmon interaction in doped epitax- 
ial graphene results in an increase in the electron scat- 
tering rate in a narrow energy region where uj ~ Ep>- 
However, at deeper binding energies a nearly linear in- 
crease of has been demonstrated with a slope of ~ 
0.025 A eV _1 , which is comparable to our measured 
value of /3. 

Because of the unique Dirac Fermion behavior and two- 
dimensionality of graphene, there has been much discus- 
sion of many-body physics that would lead to lifeti me 
broadening in ARPES measurements of graphen d 5 * 53 ^ 5 - 5 -! 
In conventional bulk crystals Fermi-liquid theory predicts 
the decay of a photohole through creation of an electron- 
hole pair to result in a lifetime which scales as (E — Ep) 2 : , 
in proportion to the number of excitation pathways that 
satisfy momentum and energy conservation. However, 
the linear dispersion of the graphene bands along with 
the vanishing density of states at Ep modify this picture. 
Hence, undoped graphene is expected to show anoma- 
lous marginal Fermi-liquid behavior, characterized by a 
lifetime that scales as (E — Ep) -1 )^ Electron-phonon 
interaction has also bee n sh own experimentally to lead 
to linewidth broadening P^J However, the interaction is 
limited by the phonon dispersion to within 140 meV of 
EpW Coulombic interactions, however can affect scatter- 
ing rates for electrons well below Ep. As noted above, 
elastic scattering due to short range correlated impuri- 
ties such as adatoms, dislocations or corrugations has 
also been shown theoretically to produce a (E — Ep) -1 
dependence on lifetime.^ 

As discussed above, prior STS measurements have con- 
firmed this linear increase for a small range of energies 
(~150 meV) in the vicinity of the Fermi level for exfoli- 
ated graphene on graphite! 4 ^ Our measurement confirms 
that this behavior persists as far as 2eV below the Fermi 
level; a log- log plot of vs uj displays a slope of ~1. 
As noted above, such marginal Fermi-liquid behavior has 
also been observed by femtosecond time-resolved photoe- 
mission spectroscopy on graphite.^ 

We now return to the topic of asymmetry in MDC peak 
shape. As many recent theoretical studies have pointed 
out, the commonly made k-independent approximation 
for Im[S(k, uj)] is not full y vali d in graphene as the dop- 
ing level approaches zero ! 53 ! 54 ! The vanishing density of 
states at uj = Ep along with graphene's linear dispersion 



near Ep places a kinematic restriction on the available 
phase space for electron-electron scattering. The scatter- 
ing pathway e~ — >• e~ + e~h + is only available for off- 
shell electrons for which uj > vp\k\ and is kinematically 
forbidden when uj < vp\k\. Thus, one expects a disconti- 
nuity in Im[S(k, uj)] at uj = k and decay due to electron- 
electron interacti on m ay be indicated by asymmetry in 
MDC peak shape! 52 * 54 ! As mentioned above and indicated 
in Fig. [7| in MDCs taken through the K point for mono- 
layer graphene additional spectral weight is present in 
the uj > vp\k\ regime. In principle, a full deconvolution 
of the ARPES intensity would recover the exact function 
form of A(k, uj). However, such a procedure would re- 
quire use of the full 2D integral specified by Eq. [2j which 
is beyond the scope of the work presented here. 



V. SUMMARY 

Photoemission on thin sheets of 2D crystals is expected 
to grow in importance as interest in single layer insula- 
tors and semiconductors increases. We have performed 
ARPES on a 2D suspended surface with well defined sur- 
face corrugation. By comparing our work with our prior 
results obtained from diffraction measurements on corru- 
gation in suspended graphene sheets 14 we have developed 
a model for understanding the effect of corrugation on 
ARPES spectra. By analyzing results obtained with dif- 
ferent photon excitation energies, we have estimated the 
contribution of surface roughness to broadening. Thus, 
despite the surface corrugation in the graphene layer, it 
is still possible to develop insights into graphene physics. 
In particular, we have shown that exfoliated suspended 
graphene is essentially undoped in its pristine form. Ad- 
ditionally, we have shown that the band structure has no 
significant deviations from linearity in the vicinity of the 
Dirac point. Our measured Fermi velocity is comparable 
to results obtained on supported graphene by transport 
and optical measurements. Finally, we have also shown 
that undoped exfoliated graphene behaves as a marginal 
Fermi-liquid with an anomolous carrier lifetime, which 
scales as (E — Ep)~ x . 
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Appendix 

In this appendix, we will follow the standard formalism 
for single photon photoemission using the dipole approx- 
imation. We will adapt the treatment to deal with a lo- 
cally curved surface using a specific initial state described 
by the tight binding model for graphene. 

According to the standard tight-binding scheme, initial 
7r states in the valence band of graphene with energy 
ujk and crystal-momentum k are represented as a linear 
combination of molecular p z orbital states: 



Mr) 



iE^ E Cf^r-R-rj) (A.l) 

ViV R j=A,B 



where -4= is an overall normalization factor, A and B 
designate the sublattice sites and Tj their locations within 
the unit cell. The Cf are complex coefficients obtained 
from the tight-binding model and the (j) are molecular p z 
orbitals. The sum over R runs over all N unit cells in the 
crystal (note that we work in the limit where TV — >• oo). 

The transition-matrix for photoexcitation from this 
initial state to a plane wave final state with total mo- 
mentum kf outside the crystal can be written, using the 
dipole approximation, as follows: 



M k oc J d 3 re-^f r (p ■ A)^ki(r) (A.2) 

Inserting the above definition for the initial state we 
obtain: 



M 



cx(k i .A)^e i(kl - kf) " R 0(fc / ) <^ k e~* kf " Tj 



R 3=A,B 

(A.3) 

where <j) is the Fourier transform of the molecular 
p z orbital and A represents the polarization vector of 
the incoming radiation. We are interested in a small 
region of momentum-space in the vicinity of the K 
point. Since ki • A and 0(kf) are nearly constant 
in this region, we concentrate our attention on the 



sum, ^2 e 

EC 1 



i(ki-k f )-R 



EC 



iki-r,- 



The sum over j, 



ki -ikf-Tj 



depends only on the relative phase be- 



tween the C k, s and the pathlength difference from atoms 
A and B to the detector. This term changes rapidly on a 
contour around the K point. Along the TK direction, the 
term changes from 2 to as we pass through the K point 
from the first to the second BZ. However, if we restrict 
ourselves to the region of k-space along the TK direction 



pact 

within the first BZ (see Fig pi, ^C ki e _zkf Tj is nearly 
constant. Thus, we are left with the sum ^ e *(ki-k f )-R^ 



For a 3D crystal with perfect translational symmetry 
R = niai + 7i2a 2 + n3a 3 where m, ^2,^3 are integers 
and ai,a2,a 3 are primitive lattice vectors. In this case, 
the sum over R reduces to momentum preserving delta 
function 5(kf — ki — G) where G is a reciprocal lattice vec- 
tor. However, since graphene is a two-dimensional lattice, 
momentum conservation does not hold in the perpendic- 
ular direction. More significantly, exfoliated graphene is 
a flexible membrane that is not atomically flat, so the 
R's must be expressed in terms of a continuous variable; 
thus R = niai + 712^2 + Ax + Ay + 2, where zj is a 
continuous variable which represents the local height of 
the graphene sheet. The variation in height is such that 
we can consider the well known theory of scattering from 
continuous rough surfaces in order to evaluate the sum 



in equation A.3 We begin by replacing the discrete sum 
with an integral: 



R 



i(ki-k f ) 



R = J d 3 rD(r)t 



,i(ki-kf)-r 



(A.4) 



where D(r) is the density- function of the material which, 
for a perfectly crystalline flat sample, has the form: 



D(r) = £<5(r-R) 



(A.5) 



R 



A periodic lattice generates a photoemission spectrum 
with the periodicity of the reciprocal lattice. However, 
since we are concerned with the photoemission spectrum 
in a small region of k-space in the first Brillouin zone, we 
may abandon the description of the surface as a discrete 
lattice and replace it with a smooth, continuous sheet. 
Thus, we approximate D(r) as a surface density function: 



D(p,z')^5[z'-z(p)} 



(A.6) 



where p = ry = (x,y) and z' = r±. Thus, the surface is 
now defined by the height function z = z{p). Inserting 
the above definition of D(r) and explicitly separating ki 
and kf into components parallel and perpendicular to the 
surface we obtain: 



M oc d 2 pe 



i(ki 



-kf_L)-*(p) e i(ki|| 



-k f 



(A.7) 



We have retained momentum conservation; for a constant 
z(p) the above integral produces 5(kiy — kf y) times a com- 
plex phase; but for a non-trivial z(p), the delta function 
broadens since e^ ki ^ -kf ^ ,2; ( p ) is no longer independent 
of p. Additionally, electronic states in graphene propa- 
gate on a curved space, which implies that the direction 
of the initial state wavevector, ki, varies as a function 
of position along the surface so that ki± and vary 
with p as well. This introduces additional phase varia- 
tion into the exponential argument (ki_L — kfj_) • z(p). 
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However, this variation is very small in comparison to 
that introduced by changes in z(p) and can effectively 
be ignored with little change in our final result. In par- 
ticular, ki± varies proportional to which is on the 

order of 0.01 A -1 . Thus the phase variation in the term 
ki±z <C 1 (Az ~ 2 A) is very small in comparison to 
the variation in the kf±z term (fc/j_ ranges from 2.5 to 
3.5 A" 1 ). Thus, we will approximate the direction of 
the initial state wavevector, ki, as constant for all points 
on the surface. This means that k is not p dependent 
and we can define a new vector k = ki — kf , so that our 
expression becomes: 

Moc J d 2 p^ z ^^wP (A.8) 

To find the photoemission intensity we use Fermi's 
golden rule which yields: 



that the full width at half maximum (FWHM) of I(k\\) 
(with k± held constant) scales as £ > ~ 1 (wk±)i when 
(wk_\_) 2 ^> 1. The functional form of I(kn) is well ap- 
proximated as!^ 

5 fcx (k,|) = ^/(wk ± )^)F a (k^/(wk ± )^) 

F a (Y) = J XdXexp(-X 2a )J (XY) (A.14) 

The above discussion began with the assumption 
that the initial state, i/ik, had a well defined pseudo- 
momentum, /c, and energy uo^. To include initial state 
broadening in our description, we replace ipk with a sum 
over multiple momentum states, J^a^k: where the a& 
are complex coefficients. The coefficients, a^, are related 
to the spectral function, A{k,uj) by \a^\ 2 = A(k,uS) with 
the spectral function defined as: 



iaf ||M|| 2 



||M|| 2 oc J d 2 p'd 2 pe i ^ z ^ 
xe ik ii •^ e _ik ^ ,2 (^ / ) e _ik ii 'p' 



(A.9a) 



Defining r = p — p' we can rearrange to obtain: 



I OC J d 2 V (^J ^ pe ik^.[ 2 (p+r)- 2 (p)]^ e ik r r ( Al0 ) 

The term inside the parenthes is the height-difference 
function, C(r, k±), of the surface which is related to the 
height-height correlation function, H(r) = (\z(ro + r) — 
z(ro)\ 2 }. It is straight forward to show that the C(r, fcj_) 
equals e~^ H ^ kl - ^ Thus, we have: 



Joe J d 2 rC(k ± ,r)e i ^-r^ 
C(fc ± ,r)=c-i H W fc i. 



(A.lla) 
(A.llb) 



For a large class of surfaces, H(r) has the following 
properties: 



H(r) oc e 2a ,for r <C £ 
H(r) = 2w 2 ,for r > £ 



(A.12) 
(A.13) 



where a is a measure of the small scale roughness termed 
the "roughness exponent." In, particular, it can be shown 



ira(E) 



#e(E)) 2 +/m(E) 2 



(A.15) 



where E = E(fc,cj) is the quasiparticle self-energy. Re- 
taining our simple description of the final state as a free- 
electron state with momentum q, our transition matrix 
becomes a sum, M = J^a^M^, over multiple matrix 
elements weighted by the complex coefficients a&, where 
the M kq are the original transition matrix elements de- 



fined in Eq. A. 3 Again, using Fermi's golden rule we 
find that the transition probability is proportional to the 
square of this sum: 



Icx\M\ 2 = \J2a k M kq \ 2 = 

k 

\a k \ 2 \M kq \ 2 + ^ aia k >M kq *M k ' q (A.16) 

k k^k' 

The cross terms have the form: 



M kq *M k ' q oc J d 2 r (y d 2 pj^ z ( p+r 

(A.17) 

where Aky = k'y - ky, k ± « k' ± . The e lA ^w p factor in 
the p integral introduces a random phase that causes the 
integral to average to zero (since it is taken over the whole 
surface). Thus, the cross terms can be safely neglected 
and we arrive at the final expression or photoemission 
intensity as a function of k\\ with k± fixed, described by 
Eq.0 
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